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Abstract 

The  steady,  zonally  symmetric  motion  of  a  shallow 
incompressible  atmosphere  on  a  rapidly  rotating  earth  subject 
to  an  equator-ward  temperature  gradient  is  studied.   The  assump- 
tions made  allow  the  thermodynamics  to  be  treated  separately 
from  the  motions o   The  turbulence  terms  are  modelled  using  mixing 
length  arguments  o   Ass\iming  the  turbulence  length  scale  is  small 
compared  to  the  earth's  radius  and  the  motion  is  slow  compared  to 
the  earth's  rotation  speed,  the  north-south  geopotential  gradient 
drives  the  eastward  winds  in  the  "interior"  temperate  regions, 
that  is,  not  near  the  surface,  equator  or  poles.   The  meridional 
winds  in  the  interior  are  driven  by  the  turbulence  generated  by 
the  shear  in  the  eastward  winds.   Near  the  equator  the  advective 
terms  become  comparable  with  the  rotational  terms,  but  the 
turbulence  terms  remain  unimportant o   The  motions  there  show 
trade  winds  at  the  equator,  changing  to  eastward  at  a  predicted 
latitude  of  about  15°.   The  meridional  motion  takes  the  form  of 
Hadley  cells  with  rising  at  the  equator  and  sinking  again  at 
about  21°.   The  Hadley  cell  and  the  temperate  region  are 
connected  through  a  vertical  layer  of  turbulence.   Near  the  poles 
the  advective  and  turbulent  terms  become  comparable  with  the 
rotation  terms.   Surface  Ekman  layers  complete  the  picture. 


Ill 


1.  Introduction 

When  viewed  from  outer  space,  the  earth's  atmosphere,  as 
evidenced  by  the  cloud  patterns,  shows  no  signs  of  symmetry  or 
regularityo   However,  analysis  of  wind  and  temperature  data  over 
long  periods  of  time  shows  that  on  the  average  the  atmosphere 
behaves  in  a  regular  manner.   Prevailing  westerlies,  trade  winds 
and  the  Hadley  circulation  are  gross  features  of  the  motion 
which  are  observed  from  average  data. 

In  1735j  George  Hadley  initiated  the  idea  that  solar 
heating  at  the  equator  forced  the  air  there  to  rise,  and  hence 
he  conjectured  that  the  air,  once  aloft,  travelled  to  the  poles, 
where  it  sank  back  to  the  surface  to  journey  toward  the  equator 
again.   Subsequent  observations  by  Ferrel  around  the  1880' s 
suggested  that  the  air  rising  at  the  equator  did  not  travel  to 
the  poles,  but  sank  back  to  the  surface  at  about  30  north  or 
south  latitude.   In  addition,  the  cold  air  sinking  at  the  poles 
rose  again  at  about  6o  north  or  south  latitude.   In  between 
these  cells  was  a  single  cell  rising  at  6o  and  sinking  at  30 
in  each  hemisphere.   The  conjectures  mechanism  to  drive  the 
temperature  cell  was  friction  between  the  adjacent  thermally 
driven  cells;  while  the  eastward  winds  therein  were  the  result 
of  the  rotational  (Coriolis)  force  interacting  with  the  meridio- 
nal motions.   Because  of  the  appearance  of  three  cells  in  each 
hemisphere,  mechanisms  to  explain  the  motions  by  explaining  the 
presence  of  each  cell  have  been  called  tricellular  theories. 

The  aforementioned  tricellular  theory  of  driving  the 
motions  has  fallen  into  disfavor,  one  reason  being  that  if  the 


heating  drives  the  motions  in  the  tropical  and  polar  cells,  the 
winds  in  the  temperate  cell  are  necessarily  weaker  since  it  is 
driven  by  the  others.   But  the  prevailing  westerlies  observed  in 
the  temperate  cell  are  much  too  strong  to  be  secondary  or 
frictionally  forced  winds. 

A  recent  popular  explanation  of  the  driving  mechanism  is 
the  wave  theory  (Pfeffer,  1964).   In  this  theory,  the  north-south 
temperature  gradient  induced  by  solar  heating  drives  wave 
cyclones  and  anticyclones.   The  nonlinear  interactions  of  these 
waves  cause  a  mean  zonal  flow.   The  zonal  flow  in  turn  drives  the 
meridional  cells  through  the  mechanism  of  the  Coriolis  force. 

We  give  a  new  explanation  for  the  general  circulation  in 
which  the  north-south  temperature  gradient  drives  the  zonal 
motions  in  the  temperate  cell  through  the  geopotential  gradient. 
The  zonal  motions  in  the  tropical  and  polar  cells  are  driven  by 
the  zonal  motion  in  the  temperate  cell  through  friction  between 
the  cells.   More  precisely,  the  eastward  motion  at  the  poleward 
edge  of  the  tropical  cell  causes  an  eastward  motion  inside  the 
tropical  cell,  at  least  for  some  distance.   The  Coriolis  effects 
cause  a  strong  meridional  motion  and  the  trade  winds.   The 
meridional  motions  in  the  temperate  cell  are  driven  by  frictional 
interactions  with  the  mean  zonal  flow  and  the  tropical  and  polar 
meridional  motions. 

In  order  to  describe  the  motions,  we  shall  use  the  fluid 
mechanical  and  thermodynamical  equations  for  the  zonally  sym- 
metric (i.e.  independent  of  longitude)  flow  of  a  shallow  layer 
of  fluid.   Moreover,  following  Saltzman  (I968),  we  shall  make  a 


Boussinesq-like  approximation,  neglecting  density  variations 
except  in  the  vertical  momentum  equation.   The  shallow  atmosphere 
assumption  allows  us  to  use  the  hydrostatic  pressure  equation 
for  vertical  momentum  balance.   This  assumption  also  allows  the 
neglect  of  the  vertical  velocity,  except  in  terms  involving 
vertical  derivatives,  which  are  large.   The  hydrostatic  pressure 
equation  allows  us  to  use  the  pressure  as  the  vertical  coordinate, 
replacing  z,  the  height  above  sea  level.   The  relevant  vertical 
"velocity"  is  aj  =  dp/dt,  the  material  derivative  of  the  pressure. 
We  shall  refer  to  cu  as  the  vertical  p-velocity,  and  to  w  =  dz/dt 
as  the  vertical  z-velocity.   The  pressure  gradient  terms  are 
expressed  in  terms  of  the  geopotential,  defined  as  0  =  gz((f),p), 

y^      ^      ^  /N 

where  z(<p,p)  is  the  height  of  the  isobar  surface  of  pressure  p. 

The  average  motions  which  we  wish  to  describe  are,  to  a 
large  degree,  both  steady  and  zonally  symmetric.   Thus  we  shall 
assume  that  time  and  longitudinal  derivatives  may  be  neglected. 

When  the  equations  are  appropriately  nondimensionalized, 
there  appear  two  important  parameters,  the  Rossby  number 
Ro  =  U/20a,  and  the  Ekraan  number  ^   =  e/20U.   Here  U  is  the 
velocity  scale,  O  is  the  rotation  rate  of  the  earth,  a  is  the 
radius  of  the  earth,  and  e  is  the  scale  of  the  turbulent  flirx  of 
momentum  per  unit  mass.   In  the  earth's  atmosphere,  the  Ekman 
number  and  the  Rossby  number  are  both  small,  typical  values  being 
^^  10"-^,  Ro  ^  10"  ,   The  presence  of  these  small  parameters  in 
the  problem  allows  the  hope  of  meaningful  approximate  solutions 
valid  in  different  parts  of  the  meridional  ( 9  -  p )  plane. 


For  example,  the  Ekman  number,  which  measures  the  impor- 
tance of  turbulence,  is  quite  small,  and  presumably  multiplies 
the  highest  derivatives.   This  implies  that  the  regions  where 
turbulence  is  important  are  boundary  layers  which  are  thin  com- 
pared to  the  scale  of  global  motions.   The  Rossby  number,  which 
measures  the  importance  of  advection,  is  also  small,  although 
not  as  small  as  the  Ekman  number.   Thus  the  advection  terms  are 
important  only  in  relatively  small  regions  compared  to  large 
scale  motions.   The  regions  where  advection  is  important,  how- 
ever, are  large  compared  to  the  turbulent  or  Ekman  layers. 


2o  Equations  of  Motion 

The  equations  of  motion  and  heat  balance  for  the  zonally 
symmetric  flow  of  a  shallow  atmosphere  (Saltzman,  I968)  are 

(2.1)  ^4-  (f  -^^^^)v  -  X  =  0, 

^P        a  cos  (f)  S(j) 

(2.2)  V  -^  +  ^  ^  +  (f  +  i^  a)G  +  ^  -  Y  =  0  , 

aS(|)     ^P  aS(j) 


(2.3)        4  +  X  T  =  0  , 
^P   P 


(2.4)        ^  +   ^^  '^"^  ^   =  0  , 


^P   a  cos  (f)  ^(j) 


(2.5)      V +  (-^ -)(Si    _  Q  =  0 

a^tj)    ^P   c  p 


•V       /\ 


A  caret  is  used  to  denote  dimensional  variables.   Here  u,  v  are 
the  eastward  and  northward  components  of  velocity,  tj)  is  the 
latitude,  p  is  the  pressure,  cd  =  dp/dt  (t  is  time),  $  is  the 
gravitational  potential,  a  is  the  planet  radius,  T  is  the  temper- 
ature, R  is  the  gas  constant,  and  c   is  the  heat  capacity  of  the 
atmosphere.   Also,  f  =  20  sin  tp  is  the  Coriolis  parameter,  where 
O  is  the  angular  velocity  of  planetary  rotation.   The  quantities 

/S       /N       /N 

X,  Y,  Q,  to  be  specified,  represent  turbulent  diffusion  and 
thermal  forcing  of  the  atmospheric  motions. 

It  is  not  our  goal  to  study  the  thermodynamics  of  the  atmo- 
sphere.  ¥e  shall  assume  that  we  can  prescribe  the  gravitational 

^   ^   /^ 

potential  $((}),p).   This  assumption,  along  with  the  Boussinesq 

approximation,  effectively  separates  the  thermodynamics  from  the 
mechanics,  since  by  determining  the  motions  in  terms  of  $,  we  can 
eliminate  the  velocity  components  from  the  thermodynamic  equa- 
tions,  yielding  equations  for  $,  T  and  other  appropriate  vari- 
ables.  Saltzman  (I968)  discusses  a  similar  approach. 

The  region  in  space  where  these  equations  govern  the  motion 
is 

-  2  -  ^1  2 

0  I  P  ^  Ps^^^  ' 

where  p  {^)    is  the  pressure  at  the  surface. 

The  boundary  conditions  which  are  needed  to  complete  the 
system  fall  into  two  categories,  turbulent  (viscous)  and  symmetry. 
The  turbulent  boundary  conditions  are 


(2.6a,b) 
(2.7) 


/s       /x 


u  =  v  =  0   on  P  =  Pg((t)), 


no  stress   on  p  =  0 


Because  the  level  p  =  0  actually  corresponds  to  z  =  oo ,  an 
unrealizable  situation,  we  shall  not  attempt  to  satisfy  (2.7). 
This  bo-undary  condition  can  presumably  be  satisfied  by  inclusion 
of  a  turbulent  boundary  layer  at  p  =  0.   A  more  realistic 
approach  to  the  boundary  condition  at  the  "top"  of  the  atmosphere 
might  be  to  apply  a  condition  at  some  p  =  p-^   «   Pg(4')  corre- 
sponding to  a  shallow  atmosphere,  say  at  about  50  km.   The 
appropriate  condition  should  reflect  the  presence  of  the  atmo- 
sphere above  50  km. 

The  symmetry  boundary  conditions  are 


(2.8a,b) 
(2.9a,b) 


V 


Su 

hi 


=  £H=o  at  ({)  =  0, 


V  =  u 


Tf 


0   at   ^  =  ±  -^ 


Let  us  nondimensionalize  the  equations  with 


u  =  Uu  , 


(2.10) 


\ 


V  =  Vv  , 


P  =  P^P  , 


(I)  =  i'^  + 'i>^'l>  . 


0^  =  (pQV/a(t)^)a3  , 


0  =  (J)  $  , 


(X,Y)  =  0(X,Y)  , 


where  U  is  a  velocity  scale,  cj)   is  the  latitudinal  scale  of  the 
motion,  p   is  a  pressure  scale,  taken  to  be  the  surface  pressure 
at  the  equator.   The  scale  of  co,  p  U/acf)  ,  is  suggested  by  the 
continuity  equation  (4).   We  also  assi;iine  that  $   is  an  appropri- 
ate scale  for  the  gravitational  potential,  and  6  is  a  scale  for 
the  turbulent  stresses. 

We  shall  denote  e/aDU  by  [   ,  and  U/20a  by  Ro.   We  call 
these  two  parameters  the  Ekman  number  and  the  Rossby  number, 
respectively.   It  is  important  to  note  that  the  Rossby  number 
defined  here  is  not  that  used  by  Charney  (1948,1963).   He  includes 
the  quantity  sin  cj)  in  his  Rossby  number.   The  Rossby  number  we 

As. 

have  defined  is  a  constant,  independent  of  ^, 

In  terms  of  these  parameters  the  nondimensionalized  equa- 
tions become 

„  /   ^u  cos((})  +  A  (j) )     -^  \ 

2 
(2.12)   HR2  (V  |^  +  cd|^)+Ro  tan(4)^  +  (^^(^)u2  +  sin((t)^+^^(f))u 


.     hY    cos(c[,  +cl)  i) 

2.13)    ^  +  n ^     ,  s^.   =  0  . 

^   cosU)  +(^^(^  ^c^ 


c   ^o 


3 •  Turbulence  Model 

Let  us  now  consider  models  for  the  turbulence  terms.   The 
northward  transport  of  eastward  momentum  due  to  turbulence  is 


(3.1)  X  = ^—j  A  [(v'u'  )  cos  J]  , 

a  cos  (j)  Stp 

where  a  prime  denotes  the  fluctuation  from  the  zonally  averaged 
flow,  and  the  bar  denotes  the  zonal  average.   We  shall  assume  a 
mixing  length  hypothesis  based  on  the  assumption  that  differences 
in  the  eastward  velocity  are  most  effective  in  supporting  the 
turbulence.   For  a  discussion  of  mixing  length  theory,  see  Hinze 
(1959). 


We 

assume 

that 

(3.2a) 

a- 

and 

(3.2b) 

V' 

^2   hn 

where  i-,  and  i      are  mixing  lengths «   The  northward  transport  of 
eastward  momentum  due  to  turbulence  is 


2  2 


bu 


b^ 


cos  ({)} 


where  L„  is  the  horizontal  size  scale  of  the  eddies,  and  Att  is 

-M.  11 

an  average  dimensionless  mixing  length,  defined  by 

(3.4)  l|a|  =    I  1^1  . 

The  appearance  of  the  minus  sign  is  accounted  for  by  the 
following  argument.   If  ^u/^(j)  >  0  and  £-.    >  0,  then  u'  >  0;  but 

the  average  momentum  transport  is  southward,  since  we  are  moving 

1  Su  ^ 

a  parcel  with  excess  momentum  -p^ —   southward,  so  that  v'  <  0. 

.   -      ^^ 

Similar  arguments  for  Su/^c})  <  0  lead  to  (3.3). 


The  vertical  transport  of  eastward  momentum  is 


(5.5)  X^  =  -  ^  (cd'uM  o 

If  we  again  use  mixing  lengths,  we  have 

(3.6a)  a-  =  ^^ 

and 

(3.6b)  ^  =  pgi^  ^  , 

^P 

where  we  now  assume  that  the  velocity  differences  which  are 
effective  in  the  vertical  momentum  transport  are  vertical  differ- 
ences, so  that  the  relevant  derivative  is  ^u/Sp.   Here  H^    and  i,u 
are  effective  pressure  differences  which  scale  the  vertical 
turbulent  eddies.   The  quantity  pg  converts  vertical  z-velocity 
to  vertical  p-velocity.   This  hypothesis  leads  to  the  following 
expression  for  the  vertical  transport  of  eastward  momentum  due  to 
turbulence: 

2     2 

(5«7)  X^  =  y-   ^  [P\(^)  ]  . 

^o 

Here  K  is  the  vertical  size  scale  of  the  eddies  and  A  is  an 
average  dimensionless  mixing  length,  p  is  the  density  at  the 
surface,  and  p  =  p  p. 

For  the  turbulent  transport  of  northward  momentum,  we  shall 
use  similar  considerations.   Again  we  shall  assume  that  differ- 
ences in  eastward  velocity  are  most  effective  in  supporting  the 


turbulence.   This  leads  to  the  expression 

2  2 

/.ON  0       V    a    r   2,^Ux2. 

^o 
and 

2  2 

^o 

where  \j,„   and  ii.„  are  average  dimensionless  mixing  lengths.   ¥e 
H       V 

shall  assume  that  'K„,    A„,  \i„   and  u^  are  constants.   Then,  without 
loss  of  generality,  we  can  take  Att  =  !• 

The  correct  scale  8  for  the  turbulence  terms  is  the  larger 
of  the  vertical  and  horizontal  scales, 

^2    ,,2   .2„2 


(3.10)  e  =  max  ) ^ — '  ^IT 

1    Po    ^  K 

By  taking  L„  =  100  km,  and  L^.  =   100  mb,  corresponding  to  a 

vertical  eddy  scale  of  a  couple  of  kilometers,  we  find  that  the 

horizontal  and  vertical  turbulence  terms  are  comparable  if  ({)  =1. 

If  h      is  smaller,  the  horizontal  terms  are  dominant.   Unless 
^o 

otherwise  specified,  we  shall  use  the  horizontal  turbulence  terms 

for  scaling  purposes.   We  denote  L  /a  by  e .   We  now  assiime  that 

2 
e  <<  1  and  Ro  <<  1.   Since   &  =  1_  . -p£  ,  the  turbulence  terms  will 

be  negligible  compared  to  the  advective  terms  as  long  as  i      »   e. 
Thus  for  latitudinal  scales  larger  than  a  typical  large  scale 
eddy,  we  can  neglect  the  turbulence  terms.   More  precisely,  the 
role  of  turbulence  has  been  relegated  to  various  boimdary  layers 
in  the  flow. 


10 


4 .  The  Interior  or  Temperate  Approximation 

In  seeking  approximate  solutions,  it  is  natural  to  assume 
(j)  =  1,  (|)  =0,  and  a  =  1  and  neglect  terms  of  order  Ro.   This 
yields  the  approximation 

(4.1)  -  sin  ^v  ^   0  , 

(4.2)  sin  (f)u  =  -  ^/S(j)  , 

op    cos  9  d(p 

Here  we  have  set  0   =  20Ua. 

Equation  (4.1)  suggests  not  that  v  =  0,  but  that  the  first 
term  in  the  expansion  of  v  is  small  compared  to  1.   An  Ekman 

layer  expansion  for  the  surface  boundary  layer  (Pedlosky,  1959) 

1/2 
suggests  that  the  next  terra  in  the  expansion  is  of  order  ^  '     . 

Pedlosky  ignores  the  inertia  terms  and  uses  a  viscous  term  to 

model  the  effects  of  turbulence.   Inclusion  of  the  inertia  terms 

does  not  change  the  result.   However,  at  this  order  in  the  inner 

equation  v  is  still  zero. 

The  appropriate  balance  in  the  inner  region  is  a  =  G  ,  so 

that 

(4.4)  -sin  (|)v  =  X  . 

From  equations  (4.1-4.4),  the  mechanism  driving  the  motion 
in  the  interior  (temperate)  region  is  as  follows:   the  geo- 
potential  drives  the  eastward  motion  through  the  Coriolis  force. 
The  meridional  motions  v  and  cd  are  driven  by  friction;  more 


11 


specifically,  the  turbulent  stress  in  the  eastward  direction 
drives  the  northward  motion.   The  vertical  motion  is  determined 
through  the  continuity  equation. 

In  order  to  determine  the  motions  further,  it  is  necessary 
to  have  expressions  for  the  geopotential  0  and  the  density  p, 
which  appears  in  the  eastward  turbulent  stress  force  X.   Thermo- 
dynamic considerations  are  important  in  each  expression. 

Since  the  velocity  scale  U  can  be  taken  to  be  the  largest 
value  assumed  by  the  magnitude  of  the  velocity  vector,  it  is 
approximately  (in  dimensional  variables) 


:4.5) 


U  =  max 


^/^cj) 


20a  sin 


If  we  take  values  of  $  from  Schutz  and  Gates  (1971)  the  velocity 
scale  derived  from  the  model  is  about  30  m/sec,  which  gives  a 
Rossby  number  of  about  0.03.   The  eastward  velocity  in  the 
northern  hemisphere,  as  calculated  from  the  data  of  Schutz  and 
gates,  is  shown  in  Figure  1. 

From  the  meridional  motions,  using  the  mixing  length 
hypothesis  of  Section  3?  we  have 


(4.6) 

1     r   1     s 

"Su 

hu 

sin  (f)  Icos  (|)  h^ 

b^ 

where  p 

2  2  3 

Vv^  Po^  . 
=  T— g —  IS  assumed  1 

Po4 

:o  I 

e 

-=*]^pi0'h 
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800  mb 


90     80     70     60     50     40     30     20         <^(o) 

Figure  1.,   Zonal  wind  versus  latitude: 

at  400  mb  predicted;  o  observed; 

at  800  mb  -.-  predicted;  x  observed. 
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5.  The  Temperate  Surface  Boundary  Layer 

There  are  several  boundary  layers  needed  to  complete  the 
picture.   It  is  well  known  that  an  Ekman  (turbulent)  boundary 
layer  is  needed  to  satisfy  the  boundary  conditions  at  the  earth's 
surface  (2.6a,b).   In  the  temperate  region,  where  the  advective 
terms  are  negligible,  the  analysis  of  Pedlosky  is  valid. 
Pedlosky,  however,  uses  Newtonian  viscosity  terms  to  model  the 
turbulence.   We  feel  that  the  effect  of  turbulence  can  be  better 
modelled  by  the  mixing  length  considerations  introduced  before. 
In  the  temperate  region  we  wish  to  balance  the  turbulence 
terms  with  the  rotation  terms.   This  leads  to  a  pressure  scale  of 
(p  gir   JL)  /5^  ^^^   ^Yie    following  equations  of  motion. 

(5-1)        ^ft^)"'  =-i"*-' 

(5-2)        4   4  t(|^)^l  =  sin  Ju  +  ^  (Ip^(J))  , 

(5.3)        "Sp   ^ — T  -  u  ,        ^, 

cos  ({)  ^tp  \ 

where  A-r^  was  taken  to  be  lo 

Since  the  surface  boundary  layer  does  not  shed  any  light 
on  the  global  circulation,  we  shall  not  discuss  it  further  here. 


6.  The  Equatorial  Cell 

The  approximation  ^      =  1  is  not  valid  near  the  equator, 
where  ^   is  small,  or  near  the  poles,  where  tan  ^   is  large.   To 
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find  solutions  valid  in  these  regions,  we  must  re-scale  the  lati- 
tude by  choosing  (|)   to  be  small.   The  result  is  an  advective 
"boundary  layer",  which  is  small  compared  to  interior  sizes,  but 
large  compared  to  the  Ekman  layers. 

To  obtain  an  approximation  valid  near  the  equator,  we  let 
(j)  =  (Ro)  ''^  ,  (f)  =0,  and  a  =  1,   The  equations  then  become 


(6.1) 


lf^li=°(^°''')' 


(6.2) 


(6.3) 


Rov— r  +  cu-T—  +  (bu  +  -— p 


Ro(v|J^  +  CD§^  -  ^v)  =  0(Ro5/2^ 


-  0(Ro^/^)  , 


^u 


Su 


1/2 
In  this  case  we  scale  the  geopotential  by  0  =  (Ro)  "^  (2Un) 

Thus,  to  lowest  order  in  Ro  we  have 


(6.4) 


Sv   ^  _  „ 


(6.5) 


^u  ,    ^u 
^     op 


=  h    , 


(6.6) 


Sv  .    Sv 

V  T^-r  +  CD  ^— 

^      dp 


=  -  (|)U  - 


^ 

^ 


(t)=0 


In  terms  of  the  equatorial  dimensionless  variables,  the  region 
on  which  we  should  consider  (6.4-6.6)  is 


(6.7) 


-00  <(f)<co,    0<p<l. 


Let  us  attempt  a  solution  to  the  equations  of  motion  (6.4 
6.6)  under  the  conditions  corresponding  to  "spring"  or  "fall". 
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that  is,  symmetric  forcing.   If  0  is  symmetric,  then  ^/^(j)     =  0, 

cp=0 

and   the    equations   become 

(6.8)  v||  +  a.|H=^v. 

(6.9)  -  li --!?  =  -*"  ' 

The  relevant  boundary  conditions  are  (2.8a, b).   A  turbulent 
surface  boundary  layer  is  needed  to  satisfy  (2.6a, b). 

The  balance  in  the  equatorial  region,  as  expressed  by  equa- 
tions (6.8-6.10)  is  the  balance  between  rotation  and  advection. 
The  driving  mechanism  (source  of  energy)  is  the  boundary  layer 
surrounding  the  equatorial  region,  with  the  important  driving 
being  through  a  turbulent  layer  between  the  equatorial  region 
and  the  interior  or  temperate  region. 

The  continuity  equation  (6.10)  implies  that  there  exists  a 
stream  function  ip   with 

(6.11)  v=||   and   <.=  -||. 

Since  we  assume  that  there  is  no  flow  through  the  top  (p  =  0)  or 
the  bottom  (p  =  1)  of  the  atmosphere,  the  lines  p  =  0  or  1  are 
streamlines. 

The  system  of  partial  differential  equations  (6.8-6.10)  is 
hyperbolic;  hence  there  exist  characteristic  curves.   The  equa- 
tions for  the  characteristics  are 
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(6,12)  |i=v. 


(6.13)  -^f  ^  ^  ^ 


(6oi4)  -ar  =  *^  ' 


(6.15)  dT  "  "'^^  ' 


where  t  is  a  variable  along  a  characteristic. 

We  observe  from  this  system  that  d^/dt  =  0  so  that  the 
characteristics  for  the  system  are  the  streamlines. 

From  the  two  momentum  equations  we  see  that  the  horizontal 

kinetic  energy  E  =  ^ is  constant  on  each  characteristic. 

¥e  shall  use  E  to  parameterize  the  characteristics.   This  implies 
that  ip   =  Tp{E). 

The  eastward  momentum  equation  (6.l4)  along  with  (6.12) 
give  u  =  -^  +  u  (E).   We  then  obtain  v  from  u  +v  =  2E,  giving 


IT 
dT  -  ^  -  -^V'^-TT'  "o^  -  "o 


(6.16)  ^  =  V  =±/2E-^-uJ^-u 


The  solutions  to  the  cj)  equation  (6.I6)  depend  on  the  sign 
of  u  o   In  order  to  get  north- south  symmetry,  we  must  choose 
u  <  0.   This  gives  v  =  0  at  (f)  =  0,  and  hence  allows  a  conver- 
gence zone  at  the  equator.   The  resulting  solution  for  the  cj) 
equation  (6.I6)  is 

^   t 
(6.17)  Mt)  =  ±(t)^(E)  dn(-|-,  k) 

where  dn  is  the  delta  amplitude  function,  and  i   is  related  to 


%   by 
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(6.18) 


^-"o 


and 


2/2E 


2E  , 


1/2 


*: 


It  remains  to  determine  the  location  of  the  streamlines, 
To  do  this  we  write 

P 
(6.19)  P 


J  dp  +p^((j))  , 


P,(*) 


for  each  ^   and 


1=MZ^=  (d^/dE)  (^^/^p^  ^ 


We  shall  assume  that  at  p  =  p  ((|')j  v  =  0.   Thus 


E 


(6.20)       p  = 


where 


d^/dE ' 


o^^^  ±J2E  -u^  -u^({>  -^ 


dE'  +p  (^) 


:6.2i) 


2E^(^)  =  (%(Eo(*))  +i^T  . 


^ 


We  now  impose  the  conditions  that  p  =  0  and  p  =  1  are 
characteristics.   We  note  that  since  tf)  =  0  lies  on  both  these 
curves,  the  value  of  E  corresponding  to  these  curves  is  given  by 


(6.22) 
This  yields 


2^1  =  (%(E^))' 


(6.23) 
and 


Po(4>)  =  1/2 
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% 


1  _  r     d^/dE' 

E  (i)  /^3"i    2     Ti   f^ 
o^^W2E  -  u^  -u^(t)  -^ 


(6.24)  ^=j  ^^         .  dE' 

o    o^   ^T 


This  integral  equation  gives  a  relation  between  ^(E)  and  u  (E). 
Thus  for  this  system,  on  the  domain  given  previously,  there  Is 
one  arbitrary  function,  which  we  may  take  to  be  u  (E). 

The  appearance  of  an  undetermined  function  in  the  solution 
is  a  common  occurrence  in  inviscid  flows.   Presumably  the  func- 
tion u  (E)  can  be  determined  by  including  dissipative  effects, 
as  in  Batchelor  ( 1956 ) .   The  condition  which  u^(E)  must  satisfy 
is 

Xdt  =  0  , 


E=const 

where  X  =  -2-  (-2H)  .   This  leads  to  the  following  functional 
differential  equation  for  u  : 


(6.25)  Au"  +  Bu'u"  +  Cu'^  +  Du'  +  F  =  0  , 

^    -^  ^  0000      o 

where   a  prime   denotes   differentiation  with  respect   to  E,    and 

2^+  A,,!  2 


A    =!5)(l)CD^dt/^'"     , 


B   =Wo?dt    , 


D   =   -  W)  4)0),  dt/^'    , 
F   =gS  <l)dt   , 
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are  functionals  of  u  .   Since  the  equation  is  second  order,  we 

o 

need  two  conditions  on  u  (E-,).   These  will  be  discussed  in 

o   1 

Section  ?• 

We  shall  not  attempt  to  solve  equation  (6.25)  here. 
Instead,  let  us  examine  the  motions  obtained  by  taking  u   con- 
stant.  One  implication  of  this  assumption  is  that  the  eastward 
velocity  is  then  independent  of  the  vertical  coordinate.   This 
assumption  is  somewhat  unreal  physically,  but  allows  us  to  solve 

the  integral  equation  (6.24).  „ 

u 
With  u  taken  to  be  a  constant,  E,  =  -^  and 


(u  .£)^ 


o_ 2_ 

2 

and  its  solution  is  easily  determined  to  be 


E  (<|))  =  p .  The  equation  itself  is  an  Abel's  equation. 


d^  _  2  1 


(6.26)  H^  -  - 

^       ^  dE     TT 


J  u  -  2E 
^  o 

Substituting  this  result  in  equation  (6.20)  for  p  gives 


1    p  /  2E  -  u^  -  u  i^  -  r}>^/4 

(6.27)       P  =  ^  ±  -  arctan   "   " . 

"^    ^         J       u^  -  2E 
V        o 

The  resulting  motions  are  shown  in  Figures  2  and  3' 

The  meridional  motions  show  Hadley  cells,  as  observed  in 

the  tropical  regions  of  the  earth's  atmosphere.   These  cells 

^  1/2 

extend  to  latitude  ^   =   ±2(Ro|u  |)  ^   in  either  hemisphere.   If 

we  assume  that  the  appropriate  value  of  u  is  that  which  gives 

the  global  maximum  value  of  u  at  the  edge  of  the  cell,  we  have 

u  =  -1,  and  we  find  that  the  cells  extend  to  about  21  north  or 
o 
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Figure  2.   Zonal  wind  versus  latitude  (both  scaled) 
in  the  equatorial  cells. 
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Figure  5c   Streamlines  for  the  meridional 
motion  in  scaled  variables « 
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south  latitude.   With  this  assumption,  the  change  from  westerly 

^  1/2 

winds  to  trade  winds  occurs  at  (j)  =  ±  /2  (Ro|u  |)  '^  ,  corre- 
sponding to  about  15°.   The  predicted  values  agree  reasonably 
well  with  observations. 

Using  considerations  similar  to  those  in  Section  5?  we  find 
that  the  surface  boundary  layer  in  the  tropical  Hadley  cell  is 
governed  by  the  equations 

(6-30)  ||-|f  =  °- 

The  vertical  length  scale  for  this  layer  is  (p^gL^-^^  (Ro )  ^    ) 
We  shall  not  discuss  this  layer  further  here. 


7.  The  Turbulent  Connecting  Layer 

From  equations  (6.17)  and  (6.18)  we  see  that  the  solutions 

in  the  Hadley  cell  region  are  defined  only  for 

^  1/2 

I  (|)  I  <  2(Ro|u  (E-,  )  I  )  ,  even  when  u  is  not  assumed  to  be  a  con- 
stant. It  is  this  failure  of  the  solution  which  indicates  that 
a  turbulent  boTondary  layer  is  needed  to  connect  the  Hadley  cell 
to  the  temperate  region.  Moreover,  since  equations  (6.8-6.10) 
have  no  forcing  term,  the  driving  mechanism  for  the  Hadley  cell 
must  be  movement  of  part  of  its  boundary.   We  claim  that  the 
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Hadley  circulation  is  driven  by  the  turbulent  layer  at  the  pole- 
ward edges  of  the  cells. 

The  correct  scaling  for  the  turbulent  layer  between  the 

1/2 
tropical  Hadley  cell  and  the  temperate  region  is  (f)   =  2(Ro) 

1/2 
and  cf)   =  £.   We  assume  that  e  _<  (Ro)  '^  ,   The  approximate  equa- 
tions become 

1/2 
to  order  (Ro)  '     . 

These  equations  are  difficult  to  solve.   We  shall  satisfy 

ourselves  by  assuming  that  u  is  continuous,  and  noting  that  if 

we  integrate  (7«l)  or  (7'2)  from  0  to  1  in  p,  and  from  -co  to  oo 

in  (|),  we  derive  the  result  that  across  the  boundary  layer^,  the 

quantity 


/  't>'  * 


0 

is  conserved.  Physically,  the  second  condition  states  that  the 
total  northward  turbulent  flux  of  eastward  momentum  through  the 
layer  is  a  constant. 

We  can  estimate  the  thickness  of  the  turbulent  layer  by 
noting  that  its  latitudinal  thickness  scale  is  e,    the  ratio  of 
the  eddy  size  to  the  earth's  radius.   If  we  assume  an  eddy  size 
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of  a  few  hundred  kilometers,  then  e  is  comparable  to  Ro,  the 
Rossby  number,  and  ^  is  about  2  o  Thus  we  might  expect  its 
influence  to  be  felt  for  about  10°  or  so. 


8.  The  Polar  Region 

The  mechanics  in  the  polar  region  is  somewhat  more  com- 
plicated than  that  of  the  tropical  regions.   If  we  attempt  a 
solution  analogous  to  the  tropical  advective  boundary  layer 
surrounded  by  turbulent  layers,  the  turbulent  layers  are  as 
large  as  the  advective  region.   This  suggests  that  the  equations 
appropriate  in  the  polar  region  entail  a  balance  between  rota- 
tion, advection  and  turbulence.   In  order  to  model  this,  we 
assTjme  that  ^     ^  ^'    ^o  ~  ^^'    ct  =  1  ^'^'^   e  =  Ro.   The  resulting 
approximate  equations  are 

.2 

lo      ^\  „  ^U   ,      ^U       _     UV       H   ^    r  1  /^U^2, 

(8.1)      V  T-r  +  1^  ^N—  =   V  -  — [-  -  — i-  -^  f  d)  (-^T-/   I   . 

2     ^ 

to  order  Ro. 

In  terms  of  the  polar  dimensionless  variables,  the  region 
under  consideration  is  0  >  ({),  0  <  p  <  p  where  p  =  p  (^)/p  (O). 

In  order  that  the  geopotential  be  smooth  at  the  pole,  we 
must  require  ^(^,  p)/^(|)  =  0. 
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The  equations  valid  in  the  polar  region  are  not  reauily 
treatable  by  classical  means.   We  can,  however,  estimate  the 
extent  of  this  polar  region,  which  we  assume  to  be  a  cell,  as  in 
the  tropical  case.   Since  ({>  =  Ro,  corresponding  to  a  few  degrees 
latitude,  the  polar  region  should  be  several,  perhaps  10  degrees 
wide.   This  is  somewhat  smaller  than  observations  indicate  for 
the  earth. 


9.  The  Energy  Budget 

A  description  of  the  motions  of  the  atmosphere  would  not  be 
complete  without  an  explanation  of  the  energy  budget.   The  asymp- 
totic approach  taken  here  suggests  that  the  following  process  is 
responsible  for  the  dissipation  of  the  potential  energy  caused 
by  solar  heating.   That  potential  energy  is  in  the  form  of  a 
north-south  geopotential  gradient.   This  geopotential  gradient 
drives  the  eastward  flow  in  the  middle  latitudes,  and  hence  is 
changed  to  the  kinetic  energy  of  the  mean  zonal  flow  there. 
Through  turbulence  generated  by  the  zonal  flow  in  the  form  of 
large  scale  eddies  and  their  related  vertical  overturnings,  the 
kinetic  energy  of  the  zonal  flow  is  changed  into  kinetic  energy 
of  the  meridional  flow. 

The  zonal  flow  in  the  middle  latitudes  also  drives  large 
scale  eddies  at  the  interface  between  the  tropical  Hadley  cell 
and  the  temperate  region.   The  large  scale  eddies  in  turn  drive 
the  motions  in  the  Hadley  cell,  where  the  zonal  motion  is  coupled 
to  the  meridional  motion  through  the  Coriolis  forces. 
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The  zonal  flow  also  drives  the  motions  in  the  polar  cells, 
which  in  general  are  weaker  than  those  in  the  tropical  cello   For 
a  schematic  representation  of  the  energy  cycle,  see  Figure  4. 


10.  Jupiter's  Atmosphere 

It  is  of  interest  to  apply  some  of  the  ideas  presented  in 
this  paper  to  the  motions  of  the  atmosphere  of  Jupiter, 
Jupiter's  atmosphere,  with  the  major  exception  of  the  Great  Red 

Spot,  appears  to  be  zonally  symmetric.   In  addition,  the  estima- 

—  P       _  p 
ted  Ekman  and  Rossby  numbers  for  Jupiter  are  10"   and  10 

respectively.   In  that  case,  the  analysis  presented  here  should 
be  able  to  predict  some  features  of  the  motion,  provided  that  the 
atmosphere  meets  the  other  requirements  herein,  such  as  the 
shallow  atmosphere  assumption. 

Since  a  knowledge  of  the  geopotential  is  required  for  pre- 
diction of  the  winds  in  the  temperate  region,  we  shall  not  concern 
ourselves  with  that  part  of  the  atmosphere.   Some  features  of  the 
motions  in  the  tropical  zones,  on  the  other  hand,  can  be  discerned 
from  the  value  of  the  Rossby  number.   In  particular,  the  motions 
show  a  Hadley  cell,  of  width  approximately  10  .   This  width 
agrees  extremely  well  with  the  width  of  the  equatorial  zone  on 
Jupiter.   The  present  model  also  predicts  trade  winds,  changing 
to  prevailing  westerlies  at  about  8  .   This  prediction  is  not 
observed;  in  fact  the  east-west  circulation  is  observed  to  be 
eastward. 
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Figure  4.   The  driving  mechanisms  for 
global  atmospheric  motions, 
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11.  Summary 

Let  us  summarize  the  asymptotic  scheme  we  have  presented 
for  approximating  the  motions  of  the  earth's  atmosphere.   We 
separated  thermodynamic  considerations  from  the  equations  of 
motion  by  assuming  that  the  temperature  field  is  knowno   This 
allowed  us  to  treat  the  fluid  dynamical  equations  as  forced  by 
the  geopotential. 

We  assumed  a  mixing  length  theory  for  the  turbulence  terms. 
By  assuming  that  the  eddy  size  is  small  compared  to  the  earth's 
radius,  we  found  that  the  effect  of  turbulence  is  confined  to 
relatively  small  boundary  layers. 

We  also  assumed  that  the  motions  we  seek  are  slow  compared 
to  the  speed  of  the  earth's  rotation  speed.   This  implies  that 
the  inertial  terms  are  negligible  except  in  layers  which  are  rela- 
tively narrow  compared  to  the  whole  atmosphere,  but  large  compared 
to  the  turbulent  layers.   These  advective  boimdary  layers  occur 
near  the  equator  and  the  poles,  and  take  the  form  of  meridional 
cells. 

The  dynamics  of  the  motion  can  be  summarized  as  follows. 
In  the  temperate  region,  the  north-south  geopotential  gradient 
drives  the  prevailing  westerlies.   Turbulent  eddies  drive  the 
meridional  motions.   The  eastward  flow  at  the  temperate  edges  of 
the  polar  and  tropical  cells  drives  the  motions  in  these  cells. 
In  the  equatorial  cell,  the  balance  between  the  Coriolis  force 
and  the  advection  causes  the  prevailing  westerlies  to  change  to 
trade  winds  near  the  equator,  and  also  creates  the  meridional 
cell  pattern. 
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There  is  a  turbulent  layer  connecting  the  equatorial  cell 
to  the  temperate  zone. 

The  mechanics  in  the  polar  region  is  complicated  "by  the 
fact  that  the  turbulent  terms,  the  Coriolis  terms  and  the  advec- 
tive  terms  are  all  of  equal  importance. 

Approximate  solutions  for  the  motion  are  not  found  in 
regions  where  turbulence  is  important.   Presumably  numerical 
techniques  are  needed  for  these  regions. 

To  have  a  complete  description  of  the  atmospheric  dynamics, 
we  need  to  be  able  to  derive  expressions  for  geopotential  ^   and 
the  density  p,  which  can  be  expressed  in  terms  of  the  temperature 
T  as  p  =  p/RT.   The  geopotential  and  temperature  are  determined 
through  the  energy  equation  and  the  vertical  momentum  equation, 
along  with  other  balance  laws  for  quantities  which  are  important 
in  thermodynamic  considerations,  such  as  water  vapor,  cloud 
cover,  etc.   In  general,  these  equations  involve  the  velocity 
components,  for  which  we  have  approximate  expressions  involving 
the  geopotential  and  the  temperature «   In  theory,  we  need  only 
substitute  the  appropriate  expressions  for  the  velocity  components 
where  they  appear  to  yield  a  thermodynamic  system  involving  only 
thermodynamic  variables,  independent  of  the  velocity  components. 
In  fact,  however,  the  expressions  for  the  velocity  components  are 
not  simple  expressions  involving  the  thermodynamic  variables,  and 
so  this  procedure  does  not  lead  to  easily  tractible  thermodynamic 
equations.   We  shall  leave  the  thermodynamic  considerations  for 
the  future. 
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An  approximate  description  of  atmospheric  dynamics  could 
be  useful  in  predicting  climate  changes  due  to  various  effects, 
for  example,  changes  in  cloud  cover,  solar  absorption,  etc. 
Most  of  these  changes  alter  the  thermodynamics  of  the  atmosphere. 
In  addition,  the  thermodynamics  is  affected  by  the  motion  of  the 
atmosphere.   What  we  have  presented  is  an  approximate  description 
of  the  motions  derived  from  given  thermodynamic  conditions.   This 
reduces  the  problem  to  consideration  of  thermodynamic  equations. 
Thus  it  is  felt  that  this  work  represents  a  step  toward  a  useful 
analytical  approach  to  global  climatology. 
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